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REIKA FUKUIZUMI AND ANDREA SACCHETTI 

Abstract. In this paper we consider a one-dimensional non-linear Schrodinger 
equation (NLSE) with a periodic potential. In the semiclassical limit we prove 
^ , that the stationary solutions of the Bose-Hubbard equation approximate the 

stationary solutions of the (NLSE). In particular, in the limit of large non- 
linearity strength the stationary solutions turn out to be localized on a single 
lattice site of the periodic potential; as a result the phase transition from 
superfluid to Mott-insulator phase for Bose-Einstein condensates in a one- 
dimensional periodic lattice is rigorously proved. 
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1. Introduction 

Bose-Einstein condensate (BEC) typically consists of a few thousands to mil- 
lions of atoms which are confined by a trapping potential and at temperature 
near the absolute zero; and BEC is described by the macroscopic wave function 
i\) = ip(x,t) whose time evolution is governed by a self-consistent mean field non- 
linear Schrodinger equation (NLSE) [25] : 

f)il> fr 2 
ih°-g =H^ + 7 |</f CT </> , H = - — A + V(x) , (1) 

with the normalization condition J t)\ 2 dx — 1, where % is the Planck constant, 
V(x) is an external potential, 7 € R is the strength of the nonlinear self-consistent 
potential and a > is the nonlinearity power. In fact, in the case of cubic 
nonlincarity, i.e. a — 1, Eq. dTJ is usually known as the Gross-Pitaevskii equation 
(GPE), 7 = Af4irh 2 a s /m, where AT is the number of atoms with mass m in the 
condensate and a s is the scattering length. 

One of the most important feature of NLSE is the spontaneous symmetry break- 
ing (SSB) effect associated to a bifurcation of stationary solutions. Symmetry 
breaking in ground states of the GPE, with a symmetric double well linear poten- 
tial, was considered in [2l El HOI [TB] . In particular, it has been shown that when the 
total number AT of particles is larger that a critical value Af c then any ground state 
is concentrated in only one of the two wells, i.e. the symmetry is broken; the value 
of Af c is explicitly given. In contrast, for AT smaller than Af c the ground state is 
bi-modal, having the symmetries of the linear Schrodinger equation. A transfer or 
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exchange of stability takes place at such critical value: for Af < Af c the symmetric 
state is stable, while for Af > Af c the symmetric state becomes unstable while the 
asymmetric state, which raises by bifurcation at Af = Af c , is stable. 

This relevant fact opens a new light on the interesting problem of BECs in a lat- 
tice and we expect that SSB effects, with the associated transition from delocalized 
to localized ground states, may occur. A prominent example of such a quantum 
phase transition is the change from the superfluid phase to the Mott insulator phase 
in a system consisting of bosonic particles in an optical lattice where the external 
potential is a periodic function [U [IT]. In fact, optical lattices basically are ar- 
rays of microscopic potentials induced by interfering laser beams and the resulting 

potential has the form V(xi, X2, ■ • ■ , Xd) — J2i=i sm2 (jJ Xe ) > wnere i s the 
wavelength of the laser light, corresponding to a lattice periods ai = |Af, I — 1, ...d 
[13) . A first theoretical model describing such a transition has been proposed in 
the case of cubic nonlinearities (i.e. a = 1) by the so-called Bose-Hubbard model 
i: 

H = - Y, J k ,jb]b k + ~Uj2Mnk-l) (2) 

k.J<EZ d : \k-j\=l k 

where b^ and h k are the annihilation and creation operators, [bjt,bt] = Sj and 

hk = bj.b/. is the number of bosonic atoms at lattice site k. The boson field 
operator ^ is written in terms of the Wannier function W\(x) associated to the 
first band of the periodic Schrodinger operator 

*0) = b k Wi(x-k-a) 

where fc = (fci, k%, ■ • • , kd) and k ■ a = (fciai, ^202, • • • , kdad)- Finally, U can be 
evaluated as follows 

U = / \W!(x)\ 4 dx (3) 

m Jn-d 

and 

J := J kij = - / Wi{x-k- a)HWx (x - j ■ a)dx (4) 

JR d 

is the hooping matrix element between neighbouring sites k and j such that \k— j\ = 
1 (actually it is independent of the indexes). Mean field arguments predict that 
as the strength of the repulsive interaction term U relative to the tunneling term 
J in the Bosc-Hubbard Hamiltonian is changed, the system reaches a critical point 
in the ratio of U /J, for which the system will undergo a quantum phase transition 
from the superfluid (SF) state to the Mott insulator (MI) state (see also [14] and 
the references therein, for one-dimensional problem see [26)). 

Our aim in this paper is to prove the validity of the results predicted by the 
Bose-Hubbard model in the semiclassical limit for any nonlinearity power a and 
in dimension d = 1, and show theoretically the SF-MI transition. More precisely, 
assuming that the potential V is regularly periodic, in the semiclassical limit we 
prove that stationary solutions of the NLSE ([TJ restricted to the first energy band 
are given by the stationary solutions of a discrete nonlinear Schrodinger equation 
equivalent to the Bose-Hubbard model (see Theorem [JJ. 
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We would remark that, in the case of cubic nonlinearities (i.e. a = 1), a math- 
ematical justification of this approximation has been proved in the special case 
of periodic potentials basically restricted to piecewise constant functions |21j ; the 
extension of such a result to periodic potentials like the optical lattice potential 
V(x) = V sm 2 {2-Kx/X) has been given by [22] under some technical assumptions 
on the band functions. In this paper, we justify the Bose-Hubbard model for any 
smooth periodic potentials, without particular assumptions, by making use of the 
semiclassical analysis for the linear problem [TJ [5] . We also prove that the ground 
state stationary solutions, consisting of orthonormal states on the first Bloch func- 
tion, are localized on a single well, in the limit of large focusing nonlinearity. As a 
result then the transition from the superfluid phase (for small nonlinearity strength) 
to the Mott-Insulator phase (for large nonlinearity strength) is proved (see Corol- 
lary HJ. 

Hereafter, we fix the dimension d = 1; for the sake of definiteness we assume the 
units such that 2m = 1; and the semiclassical parameter h is such that fi< 1. The 
notation cf> ~ ip means that lim^_j.o ^ — C f° r some C € R\ {0}; and (<f>, (p) — J <fiip; 
and let us denote V = d x and A = d xx even in dimension one; v denotes any 
generic small quantity. We will denote eventually a constant C which depends on 
parameters v, er.., by C V<(T „ to clarify the dependence on the parameters. 

The paper is organized as follows. In Section 2 we make use of the semiclassical 
analysis for the linear problem [T] [5] [T^l [T7] in order to construct a semiclassical 
approximation of the Bloch and Wannier functions with estimate of the remainder 
terms in the norm of the Banach spaces L p . Once we have obtained the semi- 
classical results in L p then we apply in Section 3, a standard argument based on 
the fixed point method and on the implicit function theorem in Banach spaces |20j 
proving that the stationary solutions of the NLSE can be approximated, in the 
norm if 1 , by the stationary solutions of a DNLS, which is the counterpart of the 
Bose-Hubbard model. 

2. Semiclassical results in dimension 1 

Here, we introduce the assumptions on the periodic potential V and we collect, 
and extend to the spaces L p , some semiclassical results on the linear operator H, 
formally defined on L 2 (K) by H = —ft 2 A + V. We recall that the dimension d is 
equal to 1. 

Basically, we assume that 

Hypothesis 1. V € C 2 (R) is a non-negative regular periodic potential with period 
a£R, that is V(x + a) = V(x). Let 



L 



1 1 

—a, —a 

2 ' 2 



be the lattice. We assume that V{x) has a unique non degenerate minimum inside 
the lattice, i.e.: there exists Xq G L such that 

V(x) > V(x ) = V mm = 0,VxeL- {x } , V"(x ) > 0. 

For simplicity, we put xq = 0, and we denote by Xj — ja the minima points for 
V(x) in the following. 

Remark 1. The assumption V"(xq) > it is only for the sake of definiteness. 
Actually, we could weaken Hyp. Q] by assuming that V {2m \x ) > and V'(x ) = 
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... = V^ m -^{x ) = Ofor some m > 1. In /act, we cowZd admit the existence of two 
(or more) minima points Xq and xo € L, too, provided that V^ m ^(xo) 7^ V^ m \xo) 
for some m > 2. 

It is well known that the spectrum of H is given by bands. Let b — 2-n/a be the 
period of the reciprocal lattice 



B = 



40 
2 ' 2 



For any fixed value of the quasimomentum variable n € B the spectral problem 

Hip = Ecp 

with quasi-periodic boundary condition 

tp(x + a,K) = e lKa ip{x, k) 

has purely discrete spectrum with eigenvalues E„(k); the functions E n {k) are named 
band functions, and the normalized (on the single lattice cell, with the L 2 norm) 
eigenvectors ip n (x, k) are named Bloch functions. Both band functions and Bloch 
functions are periodic functions with respect to k, that is E n {n + jb) = E n (n) and 
tp n (x, k + jb) = (p n (x, k) for any j € Z, and the spectrum of H is given by bands, 
i.e. 

a(H) = U„ eN [a„,/3„] 

where 

[a„,l3 n ] = {E„(k), k e B} 
is the n-th band; the open intervals (/3„,a„+i) are named gaps. 

Remark 2. In dimension one, we always have that fi n < a n +i < P n +x (this is not 
the case in dimension higher than 1 ); and for h small enough it is known that the 
first bands have its amplitude exponentially small and the gaps between the bands 
with width of order K; that is, for any fixed N > and any < v < Sq there exist 
a positive constant C := C^jv independent of h, and h = > such that 

\Pn - On| < Ce-^-^/ h and Ch < \a n+1 - [3 n \ < 

for any n = 1, 2, . . . , N and any h € (0, h); Sq is the Agmon distance between two 
adjacent sites defined in the equation ill]) below. In particular the first N gaps are 
all open and the n-th band is centered at the n—th eigenvalue X n of the single well 
operator Hq defined below in {5j) (see Theorem 4-3 by [17| J. 

Now, following Carlsson [5] let Hq be the Schrodinger operator formally defined 
on L 2 {R) by 

H = ~h 2 A + V (5) 

where V is the potential obtained by filling all the wells, but one with center xq = 
(see [B] for details): 

V{x) = V(x) + d J ( x ) > i ( ;E ) '=9{x- ja) 

jeZ-{0} 

where 9(x) is a smooth and positive real- valued function with compact support on 
a neighborhood of x n — and such that V(x)+8(x) > jp 2 , for any x € L, for some 
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small and positive p independent of h. Both V and V <E L°°, then H and Hq are 
self- adjoint operators with domain H 2 (R.). 

Since the bottom of V(x) is not degenerate, then the eigenvalue problem asso- 
ciated to the single-well trapping potential V(x) has spectrum a {Ho) with ground 
state 



Ai = 

since V m i n — 0. Furthermore, 



V"(x ) 



h+o(h 2 ) 



dist [Xi,<j(H )\{\ 1 }] > 2CH 

for some C > independent of h. The associated normalized eigenvector ipo{%) 
exponentially decreases (see, e.g. Thm. 3.2 [3]), hence ipo £ L°° n H 2 . From the 
energy identity 



fi 2 ||V^o|||2 < K 2 \\VM\l* + (V>o, Ho) = XiUo\\ 



L 2 



it follows that 



IV^olU- <CTT 1/2 . 



(6) 



Furthermore, WKB analysis says that ipo is localized in a neighborhood of xo = 
and it behaves like 

ip {x) = h-Wafa h)e- d ~ A{x ' Xo)/h , a(x; h) = a (x) + 0{K) (7) 

where cIa(x, xq) is the Agmon distance between x and the point xq defined as 

d A (x, x Q )= [ \Jv{s)ds 

J x 

By a gauge choice ipo(x) is a positive real- valued function, in particular, ao(x) is 
positive and at least of class C 2 (M): in fact ao(x) behaves like the first Hcrmitc 
polynomial in a neighborhood of £o = and like [^(a;)] -1 / 4 outside of the neigh- 
borhood. 

In order to have decay estimates in norm we introduce the following quantity 
instead of the Agmon distance dA- 



d p A (x,x ) 



sj[V(s)-p 2 ] + ds 



and we recall that ipo satisfies the following estimate (see eq. (5.1) in [B] where we 
choose r = 1) 



D d p A (x,x )/h 



L' 2 



+ h 



e d "^ x ' Xo)/h \J^ a {x) 



< C. 



(8) 



Now, let d A {x, y) be the Agmon distance between two points x and y defined as 



d A {x,v) 



It follows that for any v > we can choose p > such that 

(1 - v) [d A (x, xo) —v]< d p A (x, x Q ) < d A {x, x ) , Vie e 
Indeed, by construction it follows that 

d A {x, x ) — v < d A (x, xq) < d A (x, xq) , V.t e L. 



(9) 
(10) 
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If x (£ L let x = x + na where x £ L and n € Z (for argument's sake let us assume 
that n > and sc > xo). Then we can observe that 

d p A (x, Xo) = d A (x, x ) + nd A (xo, x + a) 

and 

d A (x,x ) = d A (x,x ) + nd A (x ,x + a) . 
Then, it follows that 

d p A (x,x$) > d A (x, xo) — v + nd A (xo, xq + a) — nv 

j / \ j t \ nd A (x ,x + a) 

= d A {x,x ) - nv - v = d A {x,x ) — — ■ — r u - v . 

d A {x ,x + a) 

If we remark that d A (x,x ) > nd A (xo, x Q + a) then the last inequality takes the 
form 

d A (x,x ) > (1 — v') [d A (x,x ) - v"\ 

where we set v' — u[d A (xo, xo + a)] _1 = 0{y) and v" = v/{\ — v') = 0{y). Finally, 
Inequality (J9j) follows by denoting v — > v := max[u', v"\. 
Let 

So = inf d A (xi,xj) and Sq = inf d A (xi,Xj) (11) 

then, by construction of the potential V(x), it turns out that 

So = d A (xi,Xj) and S^ = d A (xi,Xj), if \i-j\ = l, (12) 

and 

So < d A (xi,Xj) and Sq < d A (xi, Xj), if \i—j\>l, (13) 

in particular from an argument similar to the derivation of the first inequality of 
(JHJ) it follows that 

5*0 - v < S p < S . (14) 

Lemma 1. Let ipj(x) = tpo(x — Xj), j G Z. It follows that 

i. For some positive constant C, independent of h and of the index j, 

UAlp <Ch-^ , Vpe [2,+oo]. (15) 

ii. For any fixed v > 

\\Mk\\L><Ce-fr-*\&-*)-'V*J^k,Vpe[l,^]^ (16) 

for some positive constants C := C v and V 1 ,v" — 0(y), independent of ' H 
and of the indexes j and k. 

iii. Finally, let Xo( x ) be a function with compact support on xo — 0, i.e. 
supp (xo) = for some open and bounded set flo, then for any v > 

Wxoipjh? < C 'e-^ dA ^^-^ h , Vp G [2, +oo] , (17) 

for some constant C := C v and v — 0{v), independent of h and of the 
index j . 
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\\v>o\\is<C\\Vwo\\U\vo\\v S 5 = ^-. (19) 



Proof. Indeed, we have that HV'jIIi, 2 = 1- Then, from ([6]) and from Gagliardo- 
Nirenberg inequality, which holds true for any p € [2, +oo], we have the wanted 
estimate 

WjWl, < C||V^||i a ||^||^ 4 < Ch- S /* ,5 = ^. 

In order to prove (| 16[) for any v > let us introduce the function 

w (x) = eti-^WMA^f!) 

which satisfies the following inequality (see Proposition 3.5.3 by [T2"] ) 

\\Vw \\ L 2 + \\w \\ L 2 <Ce^ h (18) 

for some positive constant C :— C v . In particular, from this inequality and from 
the Gagliardo-Nirenberg inequality it follows that 

p_2 

2p 

Now, let Wj(x) — wo(x — Xj), for which (fT8|) holds true. We remark that 

d p A {x,Xj) + d p A (x,x k ) > d A (Xj,x k ) = \j - k\S$ > \j - k\(So - v) . 
We have that 

UjAWl? = || Wj ^e- (1 -' y) K(-^)+ d A(^^)]/' i || LP 

< e-^-^- k ^- v ^ Wj \\ L , P \\w k \\ L 2 P 

< Ce -[(l-u)\j-k\(S -u)-2u]/h 

hence, the inequality (ITC1) follows by defining (1 — v)(Sq — v) = Sq — v' and v" = 2v 
where v 1 ', v" — 0(v). 

Finally, in order to prove (|17p we remark that || Xo^j 1 1 ^ Chr 1 ^ for some 
C > 0. Furthermore, we have that 

IIXo^lU' = Wxoe-^-^'^Wih, 

< \\w j \\v\\ X0 e-< 1 - v W*™V h \\ L ~ 

< (7 e l/ / ?i g~( 1_iy )( 1_l/ )[ d ^( f2 o> :,: j) _i/ ]/' i 

since d A (x,Xj) > inf xe f2 d A (x, Xj) =: d A (Q,Q,Xj), for any x 6 fio, and since ([9]). 
From these inequalities and from the Riesz-Thorin interpolation Theorem 

\\f\\L,<\\f\\% p \\f\\ ( ^ 2)/p , V/ei 2 nL°°, Vpe[2,+oo]. 

then we have that 

Wxo^jWlp < Ce^ [ ' {1 - v){1 ~ u)[dA(x ^ no) ' v]+u]/h , Vpe [2,+oo]. 

Then the result follows by defining (1 — v) {dA{xj, f2o) — v ) = dA(xj,Qo) — v', where 
v = 0(y). ' ' □ 

If h is small enough the first band of H is not degenerate (cf. Remark [2]), i.e. 
Pi < a>2, and the restriction of H to the spectral subspace associated to the first 
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band can be described by means of an infinite matrix. More precisely, let II be the 
spectral projection of H on the first band, let F — II [L 2 (M)] , let 

vj = Ui/) j} j e z, 

then we have that 

Lemma 2. For any v > and j £ Z there exists a positive constant C > 
independent on h and on the index j such that 

Hj - Vj\\i? < Ce^ s °-^' h and ||V(Vy - Vj)\\ L 2 < Ce-^°- u ^ n . (20) 

Furthermore 

Ui - < Ce-^-^' n , Vp € [2,+oo] . (21) 

Proof. The proof simply follows by adapting the arguments by [6J . Indeed, the 
first estimate in equation (|2U)) directly follows by the result in Equation (5.2) by 
[5J. Concerning the second estimate in (|2T)|) we apply the result in Equation (5.2) 
by [6] to the gradient and the estimate (4.11) by [6j. Then ([20]) follows. Estimate 
(f2~Tjl is a consequence of (|20|) and the Gagliardo-Nirenberg inequality. The proof is 
so completed. □ 

Let 

d\{j,k) = d A (x J1 x k ) 

and 

d A {h k)= inf [d\(j, m) + d\(m, n) + d\(n, k)] 

where, by (fT2|) and (fT3|) and by construction of the potential, it turns out that 

d\(j, k) = \j - k\S and d\{j, k) > max[2, \j - k\]S . (22) 

The set {vj}j^z form a basis of F such that (see the Main Theorem in [SJ adapted 
to a regular lattice) for any v S (0, So) small enough and fixed there exists a positive 
constant C :— C v independent of the indexes j and k such that 

Kv^Vk) - <J*| < Ce- {1 - u)d ^' k)/h . (23) 

Then 

= 1 + ^4 

where A is an Hermitian infinite matrix from I 2 to £ 2 with 

||A|| < Ce-<- s *-"V n . (24) 

Let i? = {bj t k)j,k be the inverse square root of 1 + A; the elements bj,k of the 
matrix B satisfy (see Eq. (5.18) by [6j) 

b jtk = ^-^,^) + o( e -^^^)/ ft ) 

= <5 3 fc + C L-\0--"')<lA(j,k)-»''V^ , (25) 

The set of functions 

ft 
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form an orthonormal basis of F and the restriction of H to F on such a basis {uj}j 
is associated to the matrix 

<$Ai + Wj , k + r Jtk (26) 

where 

Wj, k = ~ [(ipj,rk) + (rj,ipk)] , rj(x) = r (x - Xj) , r = 6 m ip , 

and 

|rj, fc | < Ce-V-^W) . 
Remark 3. Because ipj and 8 k are real valued functions then it follows that 

Wj,k — Wk.j and Wjj = . 
Here, we state some properties of the vectors Uj. 

Let T be the translation operator (Tf) (x) = f(x + a), where a is the period 
of V. Then, by construction it follows that ipj = T^tp and, since [H, T] = 0, 
[II, T] = 0, hence Vj = T^v . 

Lemma 3. Let T be the translation operator (Tf) (x) = f(x + a), where a is the 
period ofV, then Uj = T^uq. 

Proof. First of all we observe that the matrix A is an Hermitian matrix of 
Tocplitz type (i.e. = aj-k,o), indeed the elements of the matrix A arc such 
that 

$■ + cij.k = (vj,v k ) = (v ,T k - j v ) , 

that is, a,j t k = aj-k,o '■= o-j-k f° r some sequence {ae}e£z such that = ag. Hence, 
the matrix B = [1 + A] -1 / 2 is defined as the convergent series (indeed ||^4|| < 1) 

B-TcA- c - ( - 1} " (2n)! 1 ■ 

n— v ' v 

where, by means of a straightforward calculation, it follows that B(l + A)B = 1. 
Therefore, the matrix B is of Toeplitz type, since the power A n , n e N, of the 
Toeplitz matrix A is still a Toeplitz matrix. That is the elements of the matrix 
B are such that bj^ = bj-k for some sequence {bi}g^z such that b-e = bg. From 
these facts it follows that 

u = ^2 b - kVk and U J = ^2 b i-kV k ; 

k k 

in particular, it follows that 



□ 



Lemma 4. For any v > and j £ Z there exists a positive constant C > 
independent on h and on the index j such that 

|| - Vuj|| La < Ce-^ s °-^' h (27) 
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and 

Ui - u 3 \\ Lv < Ce- (So ' u)/h , Vp £ [2, +00] . (28) 

Proof. Recalling that ipj = T^'Vo, Vj = T^v and Uj = T^u , where T is the 
translation operator defined in Remark [3l then we can restrict ourselves to j = 0. 
From ([201) and (EH) then 

Mo-UoWv <Ce-<- s °-^ h 
for some constant C independent. The same estimate follows for the gradient, too: 

l|VW>o - uo)IU" < l|V(^ - v )\\ L 2 + ||V(uo - v )\\ L 2 
where it has already been seen that 

||V(^ - ^o)|k- < Ce-^-^^. 
The second term is written as 

V(w - vp) = VfcVffc - v = E [ b o,k - $0] v k 

k k 

and then 

iiv(uo-«D)ii£ a < EiVfc-^rnv^n 2 

k 

< El^-^ol 2 [2||V^ fe || 2 + 2||V(V fc -^)|| 2 ] . 
k 

Here the right hand side may be estimated as follows. 

||V^.|| L2 < cn- 1 ' 2 , ||V(«* - ^ fe )|U= < ce-^>-^/ h 

and 

E IVfc - *of < E K fc " *o| = II 1 - B \\ < Ce-( s °-"V n . 

k k 

These inequalities imply ([27]) and from the Gagliardo-Nirenberg inequality then 
(|2"51) follows. The proof is so completed. □ 

Lemma 5. The functions Uj and Vj can be chosen to be real-valued. 

Proof. We start proving that Vj are real-valued, because the functions ipj are 
positive. Indeed Vj = Hipj where II is the projection operator on the first band. 
From the Bloch decomposition formula (see eq. (2.1.22) by [20) it follows that any 
vector rp £ L 2 can be written as 

*w -£/«.«*«>*.<■)*. 

where B is the Brillouin zone, k is the quasimomentum (or also crystal momentum) 
variable and (p n (x,K) are the Bloch functions. The function a n (n) is called the 
crystal momentum representation of the wave vector associated to if) and it is defined 
as 

On(«) = / <fn(x,K)^(x)dx 
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Then, the restriction to the first band is simply given by for any vector ip G L 2 

(Hip) (x) = j <fi(x, n)ai{n)dn 
and it is real valued; indeed, 



(ILip)(x) — / <fi(x, K)ai(n)dK 



B 



y>i(x, —K)ai{K)dK 



B 



B 



(fi(x, K)ai(—n)dK 



since (see eq. (2.1.17) by 20)) ip n (x,—K) = ip n (x,K). Moreover, if ip(x) is a 
real-valued function, 



ai(— k) = / ipi(x, — K)ij)(x)dx = / ipi(x, K)ip(x)dx = I tfi(x, n)ip(x)dx — a\(n) 
Jr Jr jr 

We can apply this argument for if) — ipj. Once we have proved that Vj are real- 
valued functions then immediately follow that the elements of the matrices A and 
B are real-valued, and then, by construction, Uj are real valued functions too. □ 



Lemma 6. The wavefunctions Uj are such that: 

i. || Wj-Mfc < Ce~^ Sa ~ v ">\' J ~ k \~ lJ "}l h , j 7^ k, for any v 1 ', v" > and for some 

positive constant C > independent on the indexes j and k; 
ii- II Ej < Ch' 1 / 2 for some C > 0. 

Proof. In order to prove property i. we recall the following results by [B] (see, 
respectively, the proofs of Lemmata 5.1 and 5.2): 



and 



and 



where 



[««i(*i,0+*i(*i»0]/»ty (.)^, fc (.) 



< C 



L 1 



< c 



L 2 



< c 



(29) 



dA (xj,x) = inf [d p Axj,xi) + d p A (xe,x)] > d p A (x j} x) 
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and C is a positive constant independent of the indexes j and k. Hence, it follows 
that 



L 1 



< 



L 1 



\dA P {x 3 ,)+d A p (x k ,-)]/h^_^ k ^ 

^ AP{x >-- )+d ~ AP{Xk ^ /tl lM-)-M-M(-) 

[^(* 3 v)W(**,-)]/^.(.) _ „ (.)] [M .) _ Ufe( .)] 



L 1 



< c 



from which it follows that \\v 3 v k \\ L i < Ce-KSo-^OIj-fcl-""]/?* for any v' , v" > 
and for some C > 0. In order to get the estimate on Uj we recall the estimate 
(|25p , from this fact and by means of a straightforward calculation then estimate i. 
follows. 

In order to prove property ii. we remark that the sum y\ |tij(a;)|, if convergent, 
defines a periodic function with period a since Uj(x) = uq(x + aj); hence, we may 
restrict the L°° estimate of such a sum on the single interval with length a. Now, 
let n be any integer number and we set /„ = [{n — l/2)a, (n + l/2)a], where a is 
the period of the potential V{x), and we are going to estimate uo(x) in the interval 
/„. Then we observe that 



||w |U°°(/„) ^ ^2 \ b 0,r 



m\\Vm\\L°°(I n ) 



where the coefficients &o,m of the matrix B satisfy (|25p and where H^mll^^) can 
be estimated by means of the Gagliardo-Nirenberg inequality: 



1 1 Vrn Hl<~(/„) < C [|| Vu m || L 2 (J n ) + IKn|U 2 (/„)] 



where 



|Vu. 



m||L 2 (/„) 



rn\\L 2 



< Ch- 1 ' 2 



and where 1 1 ^rrx 1 1 i 2 ( z„ ) can be estimated as follows 

-9m (x) 



lL 2 (/„) 



e 9 ^ x) v r: 



< 



\e-*»W\\ L aa iIn) 



L 2 (I„ 



< (J e -{S -v)\n-m\/H 



where g m (x) — d,A P (x m , x)/H and where 





< 






< c 




L 2 (/„) 




L 2 (I 


) 



since ([5]) and (|2T)]) . Collecting all these results then it follows that 

H«illL-(i B )<Ce-( So -")l n -'l/ n , j?n, 
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from which property ii. follows: 



E 



< IK|U°°(I ) + II E l%'IIU°°(Jo 



< \\uo\\l-<r) + E Ce-( s °-^l^/ fi < Crr 1 ' 2 + Ce-( s °-"V n . 



□ 



Remark 4. TTie wave-vectors Uj(x) construct the first Block function tpx(x,n). 
More precisely, let Xa(x) be a function with compact support contained in an open 
set M C B XQ (So) such that < xo < 1 and it is exactly 1 on B Xo (Sq — a) for some 
a > fixed. Let \j ( x ) = Xo (x — Xj ) . Then, the vector 

01 (x, k) = E eiK ' S "Xj(*)uj(x), 

well approximates the Bloch function ipi(x,n) in the sense of Lemma 3.3 by |17) . 
Hence, the Wannier function W\ (x) defined as 

Wi(x) = y ^i(x, /c)d« 

is weZZ approximated by 

f <p l {x,K)dK = xv{x)u (x)=Mx)+O{e^ s °- a V h ) 

PI 

because of Lemma [J} 



Furthermore, the following results concerning the elements Wj.fe hold true. 
Lemma 7. Lei 

/3 = i/ |j - fc| = 1 

i/ien f3 > is independent of the indexes j and k and for any v > i/ien 

Jle-CSo+iO/fi < p < Ce-( So -"V h 

c 

for some positive constant C := C v > 0. 



(30) 



Proof. By the symmetry of the periodic potential then immediately follows that 
Wij is independent of the indexes when \i — j\ = 1, then we restrict ourselves to the 
cases of j = and fc = 1. Moreover, since ipo and d\ are positive functions then j3 
is a positive quantity. In order to get the inequality pop we have to estimate the 
quantity 



0>o,ri) = / ipo(x)ri(x)dx 

JR 

= E / ' t Po{x)6 m (x)tpi(x)dx 

ipo{x)0o(x)ipo(x - xi)dx + V] i/j (x)6 m (x)ip 1 (x)da 



m^0,l ' 
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Now, let xo = \/@o where the support of xo is contained in an open set fio such 
that d,A{xj, f2o) > §<Sb f° r any j ^ 0, then by (|TT)) it follows that 

^ / ipo(x)O m (x)i/jx(x)dx = l ^_ m (x)0 o (x)^i- m (x)dx 

m^0,l^ R m#0,l 

= ^2 Xo(x)ii- m (x)xo(x)ipi-m(x)dx 

m^O.l ^ R 

< X! IUo(-)^-m(-)IU 2 ||Xo(-)^l-m(-)IU 2 

7715*0,1 

< ( -r e 2i>/7i e -[c; A (a:_ m ,Oo)+<iA(2:i- m ,Oo)]/7!. 

7715*0,1 



for some /v < ^Sq provided v is small enough (here v = 0(v)) . Finally, from (fT6|) 
it follows that 

hence the r.h.s. of (|30l) follows. In order to prove the l.h.s. of ([50)1 let 

B X0 (l)={xeM : d A (x,x )<l} 
be the interval with center xo = and radius Z > 0, and let l\ > be such that 
B Xo {h) C supp(0 o ) and 



tl)o(x)ipi(x)Qo{x)dx > C / ^o(2 ; )V'o(x — x{)dx 

J B XQ (ii) 

In order to estimate the integral on Bo(h), let a > be fixed and let 

r Q = {a; : <2^(x, x ) + d p A (x, x\) < S + a} , 

where we underline that this set is not empty provided p, which enters in the 
definition of the function 8(x), is small enough. Then, from the WKB expansion 
of w (x; h) = h- 1 / 4 [a (x) + 0(H)] on the set [-2L, 2L] we have that 

/ t() (x)if>o(x - xi)dx = 

w (x; h)w (x - xi; h)e- {d A {x > Xo)+d * {x > Xl))/h dx 

B xo {h) 

w (x; h)w (x - xi; h)e- (d ^ {x ' Xo)+d ^ x ^ l))/h dx + 

B xo (h)nr a 

w {x; h)w {x - x i; ft) e -(^(^o)+^(a^i))A da; 

Sx ((i)\r Q 

> / w (x;h)w (x - xi;h)e-^ x ' X0 ^ +d ^ x ' Xl ^/ H dx 

JB {h)nr a 

>2c\B X0 (i 1 )nr a \h- 1 / 2 e-( s °+ a )/ h 

where the measure of the set B Xo (li) n T a is not zero and where 
c := min ao(x)ao(x — xi) > . 
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The Lemma is proved. □ 

Lemma 8. If \j—k\ > 1 then for any v > there exists a positive constant C := C v 
independent of the indexes j and k such that 

K fe | <Ce-^- k ^-^-i s °V h . 

Proof. The proof would make use of the same arguments of Lemma [7] In order 
to estimate the term (ipj,rk) we observe that is equal to the term (ipj^k^o) and 
thus we can restrict ourselves to terms of the type (?/>fc,fo), where |fc| > 1. We 
have that 

\(lpk,0m.1po)\ = \(^k- m ,Oo1p- m )\ < UxoV'fc-rnllL 2 llXoV^-mlU 2 
< (J e 2 v/ h e -[ d A{xk- m & )+d A (x_ m ,O, a )]/h 

where v = 0(y) and 

d A (x- m ,il Q ) > ^|m| - ^ S and d A (x k ^ m , Sl ) > (\k - m\ - S ■ 

Thus, we have to consider the sum 

e -[M+l*-"»l-§] So/ft ^ | fc | > ^ 

which can be estimated by Ce"[ |fc|(s ' " l/) ~3' s o]/'i^ The Lemma is proved. □ 

Remark 5. The matrix obtained collecting the elements (wj t k), for j,k £ Z such 
that \ j — fc| > 1, and the elements rj^k, for j, fceZ, defines a linear operator 

D : £ p (Z) £ P (Z) . 

This operator is bounded for any p £ [1, +oo] with bound 

\\D\\c(t>-*i>) < Ce-^ +a ^ h , 
for some positive constant § < a < So independent ofh andp, and for some positive 
constant C dependent on p and independent of h. 

3. Derivation of the Bose-Hubbard model 
Here, we consider the nonlinear Schrodinger equations 

ihd t ip = HiP + -/\iP\ 2 ' 7 ^ (31) 

where 7 € M, 

H=-h 2 A + V, (32) 

is the linear Hamiltonian with a periodic potential V{x), and | - 0| 2cr '0 is a nonlin- 
ear perturbation; 7 denotes the strength of the nonlinear perturbation and it is a 
real number, here we consider both attractive/focusing (when 7 < 0) and repul- 
sive/defocusing (when 7 > 0) cases. 

We introduce the effective nonlinearity parameter defined as 

Col 

,:= — 

where Co = ||uo 11^2^+2 ■ We remind that the set {ttjljez form an orthonormal basis 
of the space projected on the first band. 
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Remark 6. We point out that, in the case of cubic nonlinearity corresponding to 
a = 1, the effective nonlinearity parameter rj coincides, up to corrections of order H, 
with the ratio U / J between the two parameters defined in the Bose-Hubbard model 
Indeed, in the case of a = 1 we have that 

c = \\u \\l%tl 2 ~ = J IM^)I 4 ^~ / W^dx 

since Remark [^} Hence, the leading term of 7C0 coincides with the value of the 
parameter U defined by f3|). Furthermore, the dominant term of the hopping matrix 
element J defined by coincides with (3 

J = - J Wi{x - xi)HWi(x)dx = -Jwi(x-xx) [-h 2 A + V] W x {x)dx 
4>o(x - xi) [-ft 2 A + V] %j) (x)dx 



tpo(x - Xl) 



ipo(x)dx — /3 



in the limit of H small enough because of 

Hypothesis 2. [Multi-scale limit] We consider the simultaneous limit of small h 
and small 7 such that the effective nonlinearity parameter rj goes to a finite value. 

Remark 7. Since \\ipo\\^^+ a < Ch~ a / 2 and the hopping matrix parameter between 
neighbouring wells (3 is exponentially small as h goes to zero, then Hypothesis [H 
implies that 7 exponentially goes to zero as h goes to zero. 

Since V £ L°° then the Cauchy problem (|3Tj) is locally well-posed [7] and the 
conservation of the norm of ijj(x, t) and of the energy 

£{rP) = {^H^) + ^—U\\%il, 
a + 1 

holds true. Furthermore, since 7 is small enough blow-up cannot occur and, in 
fact, the Cauchy problem ([31]) is globally well-posed if I7I is small enough. Finally, 
we have a priori estimate 

\\Vip(-,t)h* < Ch- 1 / 2 and < CtT^, Vi>0, Vpe[2,oo] (33) 

for some positive constants C (see Theorem [2] and the associated Remark in [25]). 

We set ip(x,t) = e~ lXt / h <f){x), and we are interested in the associated stationary 
problem: 

A0 = #0 + 7 |0| 2ff 0. (34) 

Now, we are going to show that there exists a solution of this stationary problem 
(f> e i? 1 (]R), which can be approximated by means of a Bose-Hubbard model, that 
is by means of a discrete nonlinear Schrodinger equation (DNLS) of the form 

EF k - F j +r ] \F k \ 2 °F k = 0, ||F||^=1, (35) 

j€Z : \j-k\ = l 

where F = (i^)jez € £ 2 (Z) and 

E = ±f±. (36) 
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Remark 8. The existence of a solution of {3J$ has been investigated by many 
authors (see |20] and references therein). For example, the existence of a non 
trivial solution for the case where 7 < and A < inf o~(H) (in semi-infinite gap) is 
proved in |20] and the case where 7 e E, A > inf o~(H) and A ^ o~(H) (in a finite gap) 
was treated in |18j . No nontrivial solution exists for 7 > and A < inf u(H). The 
approximation that we will show ensures also the existence of a solution (j) G iJ x (K) 
°f j34\ l- We remark that by the standard bootstrap argument the solution <f> can 
gain the regularity, e.g., <f> G H 2 (M.). Thus, <fi vanishes at infinity and is also 
real-valued modulo the gauge transformation in case of d = 1 (see Lemma 3.7 of 
|20) ). Since we consider the problem in one dimension, it suffices to consider a 
real-valued solution. 

Remark 9. As for h35\), Weinstein proved in [27 that in the case of attractive non- 
linearities, i.e. for r) < 0, a minimizer for the associated l 2 -constraint variational 
problem, in any dimension d, 

I, = inf { - (^F.F)^) - |F,r+ 2 : £ \F 3 \ 2 = \ V \^} 



exists provided that a < -j, where 

(S 2 F) k = F i~ 2dF k ■ 

]ei d ■. b'-fc|=i 

If a > 2 an d n < 0, a minimizer exists only when \r]\ is larger than a threshold 
value ijthresh > 0. In other words, r\thresh is a value which ensures — oo < I n < 
for all r\ < with \rj\ > i]thresh- This threshold value is given by the best constant 
in the discrete Gagliardo-Nirenberg inequality, 

Vthresh = {0- + l)X, 1 = mf . 

This result implies the existence of a family of ground states of &35\) for some E. 
It is known that for E < —2d there is no nontrivial solution. 

Remark 10. In the case of d = 1 the solution F = (Fj)j e % £ £ 2 (Z) of equation i f?5)) 
can be assumed to be real-valued (see, e.g. Lemma 3.11 by [20 ); and furthermore 
it decays exponentially at infinity: i. e. there exist C > and r > such that 

\Fj\<Ce- T ^, jeZ, 

(see Theorem 1.1 in [IH]j- In particular F £ ^i(Z) C ^{Z), where 

£ 2 {Z) = J c = ( Cj ) j& , 5^(1 + |jf)| Cj | 2 < oo 

[ 3& 

Also, the case d — 1 allows us to have much more information than the case of 
d > 2. The existence of two families of positive localized modes, known as the 
site- symmetric solution, the bond- symmetric solution, is established in [24j for the 
case E > 2, with the use of a method in the dynamical system. Analysis of the 
anticontinuum limit has been done by a simple application of Implicit Function 
Theorem in [201 Theorem 3.8]; for \"q\ large enough then i35\) admits a family of 
ground state with energy E = 2— 77+0(1) in the limit of large \rj\ and with wavevector 
(Fj)~Si for any ja eZ. 
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We will need the following assumption too. Recall that we fix the dimension 
d = l. 

Hypothesis 3. Let F = (-Pj)jez be a real-valued solution of DNLS [35\) . Then, 
we assume that the linearized map at F 

L+:4(Z)^4(Z), 4:={c={ Ci }€^(Z) : Cj e K} , 

defined as 

(L + c)j = -( Cj - +1 + c J -_ 1 ) + (^ + ?7 (2 C 7+l)|F J -| 2CT )c i , c€4(Z), (37) 

is one-to-one and onto. 

The main theorem in this section is the following. 

Theorem 1. Let Hypotheses 1-3 be satisfied. Assume h > small enough and 
a > 1/2. Let E be such that \35\) admits a solution, and fix F € ^(Z) one real- 
valued solution of DNLS H35\) associated to E. Then, there exists Cn > such that 
Eg. \3J$ with \ = \\ — f3E admits a unique real-valued solution <ft G H 1 (M) close 



i.e. 



< C n , (38) 



where the constant C% — O (e a ^ h ) , for some fixed a € (0, So/2), is exponentially 
small when h — > 0. 

Remark 11. The solution <fi given in Theorem\J\is not normalized to 1. However, 
recalling that F is normalized to 1 then it follows that ||0||l 2 = 1 + C%, where 
C h = 0(e- a l h ), for some < a < So, is exponentially small as h goes to zero. 
Hence, Theorem [7] yields to a normalized solution provided that we replace 7 by 
l(l + C h ) 2 °. 

Remark 12. A similar result has been previously obtained [21] in the case of a = 1. 
In particular, in |21j the estimate of the remainder terms i38\) and J^g| ) are given 
in the energy norm defined as ( assuming that V > 0) 

\Mw ~W\\ 2 m + tl- 2 (^Vi>). 

In order to get such estimates an assumption on the periodic potential V is required; 
in fact, the results by 21 do not apply for any periodic potential, but only for 
periodic potentials with degenerate minima with infinite order (i.e. V^ n \xj) = 
for any n € N); that is practically piecewise constant periodic potentials. With more 
details, [3T] considers the NLSE in the form iut = —u xx + H~ 2 Vu + \u\ 2 u where 
the stationary solutions of the DNLS give, by means of the Wannier functions, the 
solutions of the NLSE. In order to get this result several assumptions on the band 
functions are given, in particular it is reguired that the mean value of the first band 
function of the Bloch operator —d xx + h~ 2 V is bounded when % — > 0. In fact, WKB 
arguments imply that this condition is fulfilled only when V has a minimum point 
infinitely many degenerate, if not the mean value of the first band goes as hr a for 
some a > and then the assumption above is not satisfied. In [22] the extension 
to smooth periodic potentials is given (again in the case of cubic nonlinearity) , 
provided that the band functions satisfy to a technical assumption concerning the 
cubic terms (see Assumption 2-(iii) |22j ). 



DERIVATION OF THE BOSE-HUBBARD MODEL IN THE MULTISCALE LIMIT 19 



(41) 



As a result of Theorem [1] Lemma [TU] (below) and Remark \TU\ then the phase 
transition from delocalized states, for small \r/\ (superfluidity phase), to localized 
states, for large \t)\ (Mott insulator phase) follows: 

Corollary 1. Let Hypotheses 1 and 2 be satisfied, and let a > 1/2. If the absolute 
value of the effective nonlinearity parameter n is large enough then equation \34\ ) 
has a family of stationary solutions, where each solution is localized on one single 
well, with energy A = Ai — j3E = X± — (3{2 — r\ + o(l)) in the limit of large \rj\. 

3.1. Proof of Theorem Q3 By Remark [TT)I a solution of (|34|) -if exists-is anyway 
real- valued up to gauge choices. Thus, in order to prove Theorem [TJ we look for a 
real-valued solution <j){x) of the stationary equation under the following form 

4> = + ^ , (39) 

where tfi 1 — 110 is the projection of the first band and where cj) 1 - = IT x 0, H 1 - = 1 — IT, 
is the projection of <p on the other bands. By construction, it turns out that 

a(HU) = [auPi] and a{HH^) C [a 2 , +oo) 

where 

Ch < a 2 ~ fa < ~fi (40) 

for some positive constant C . By substituting (|3"9")l into and projecting on the 
first band and on the other bands we obtain 

A0 1 = H4) 1 +7n|0| 2CT 
A0- 1 = Hcj) 1 - + 7 n- L |0| 2 ' T (/) ■ 

If we set 

then the first equation of (|4"Tj) reduces to the following system 

\ Cj = (uj,H(p x } +-/Lj, Lj := (u 3 , n|0| 2 » (42) 

where 

(uj^Hcf) 1 } = ^2c k {uj,Hu k ) . 

k 

We underline that 

ii^iU" = \wy 

because the set of vectors {uj} is an orthonormal basis of nX 2 (R). The matrix 
with elements (uj, Huk) can be written as 

{{ Uj ,Hu k )) = Xi.l-pT+D 

because of (f2"o]l , Lemmata [7] and [5] and Remarks [3] and [SJ T is the tridiagonal 
Toeplitz matrix 

m _ / o */ \j -k\?i 
(lh ' k ~ \ l if\j-k\ = l 

and the remainder term D, obtained collecting the elements Wj t k, when \j — k\ > 1, 
and rj t k, is such that 

PIU(^) < Ce-^+ a ^ n , p e [1, +oo] , 
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for some < a < So and C = C p . 
First, we justify the existence of 

Theorem 2. Assume a > and fix any number E £ R and So > 0. For any 

c = (cj)j e z, £ ^(Z) wii/i ||c||^i(g) < 6q, there exists a unique smooth map <j> : 
£ 1 (Z) x K — > i? 1 (R) suc/i i/ia£ 0" 1 = <j>(c,K) is a solution of the second equation of 
l\ l for small h > 0. Moreover, it is exponentially small as H — > in the sense that 
for any < v < So there exists a positive constant C > 0, independent of h, such 
that 

H^llffi <Ce-( s °- I/ )/ ft ||c|||f +1 . (43) 

Proof. Note that the operator H — (Ai + f3E) on IT- L L 2 has the inverse for H 
sufficiently small thanks to (|40|) . Precisely, by the functional calculus, since V is a 
bounded potential then there exists a constant C\ > independent of such that 

- (Ai + PE))- X ^\\c{l^hi) < Cih- 1 . 
Then the second equation of (|4Tj) may be written as 

J -=F(^ L ), (44) 

where we set Ai = A + and 

F(0 X ) = - 7 (tf - (Aa + pE))- l U^\<t>\ 2 °<t>, <$> = cf 1 + <£ x . 

We wish to make F be a contraction mapping in a complete metric space. With 
this aim, we make here several remarks about the property of X^ez c j u ji x )- 
First, it follows from Lemma 0] that 



\i>i 



for p € [2, oo]. Furthermore, Lemma [T] gives us the following estimate. 

\\uj\\ L v < Chr 2 ^ +Ce- {s °-^/ h . 

Remind that C does not depend on h nor on the index j. Then, for any p 6 [2, oo] 



H C J U 3 
3& 



L" 



< 
< 

< 
< 



$3i Cj -i{cfi-^ 



J 6 



We set 



,4 7 2 2ff + 1 Ci. 

Fix h > so small that (we recall that 7 is exponentially small because of Hyp. ^ 



7 2 2 "+ 1 Ci 

n 



< 1/4, 



C, II 



< 5(h), 



7 2 2<T C 



£2(2,7 + 1) 



i ll5_>H i ||l'L <l/4. 

iez 
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We shall show that F is a contraction map in 

K = {e e ^ X (M) n U X L 2 (R) : < 5(h)}- 

Indeed, for any £ = £ x + rj = -q 1 + rj X , 



11^)11^ 



< 



< 



< 



-|||sf^ll 



\\F(^)-F(^)\\ H1 < 



7^; 

7 2 2g+ 1 C; 

h 

7 2 2 "+ 1 C 1 
ft 

7d2 2 °- 



L 2 



l/||tl||2<r+l 

VIIS llj,2(2cr+l) 



|t-L||2<T+l 

IS II £,2(2(7 + 1) , 



+ -5(h) < 5(h), 



(\m%« + Millie -V^hi 



< + * 

V 4 ft 4 7 2 2ff + 1 Ci 



II? 1 - V 



IH 1 



Then there exists a unique solution ^ = <fi(c, K) € if of Eq. (j44|) for small ft > 0. 
Moreover, by the construction of the solution (see [9]), 



\H 1 



< 



for any v g (0, So) because of Remark [71 Second inequality holds true because 

III^T^Ik^ll^ll^i + llvd^l 2 ^ 1 )!!^; 

indeed, we may estimate the second term as follows. 

Iivaan 1 )^ < anv^Mie 1 !! 2 ^ 

< <? a n-w\\*\»-\\?\% 

since 



ive 1 ! 



L 2 



where we have used (|2"T|) and ^ for the second inequality, and the estimates 



1 1 1 2(7 



2rr 
£1 ■ 



for some positive constant C depending on a, complete the proof. 



□ 



The next lemma is concerned with the first equation of (1411) . The term Lj in 
(|42| may be expressed as follows 

^ = Wu^+Uc^c, + f h Si = /i(c, X ) • (45) 

We set f(c, (f)- 1 ) — (fj(c, ))j- £Z i where <f> — <f>(c, h) is given in Theorem^ 

Lemma 9. Let h sufficiently small such that Theorem^ holds. Assume a > 1/2, 
and fix any number 5q > 0. For any c = (cj)j^z G I (Z) with ||c||^i( Z ) < do, we 
have 



11%,^)^ < Cmax{4 2CT+1)2 ,C +1 }e-( s — 



(46) 
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for any v S (0, So). The positive constant C is independent of h. 

Proof of Lemma [H As a first step we remark that cf) 1 satisfies the following 



estimates 



and 



II^IU' = ||c||* <S 



(47) 



The first inequality simply comes from the inequality ||c||^2 < ||c||fi < So. Con- 
cerning inequality (|47fl we simply remark that 



liv^l 



L- 



<^|c j |||v % |U 2 <a^- 1 /2|| c | 



e 1 



L' 2 



Hence, 



Now, in order to prove the Lemma we set, for j G Z, 
/i(c,0 x ) = / i) i(c)+/ J - 2 (c,0 ± ) 

where 



/ J - 1 (c) = ( % -,|^| 2 V>-||^ll^ + 2 2 |c J f CT c j 



and 



/ J - 2 (c,^) = <u J -,|^-|^| 2 ^ 1 }. 
In order to estimate the £ -norm of (fj,2(c, <^))j g z we remark that 



(48) 



(49) 



< 



Ei 



IM a< ^-l*TVlu* 



where the first term is estimated by Lemma|S]-ii. and where, for what concerns the 
r.h.s., we remark that 

- m 2 ^ 1 = i^ 1 + ^i 2 -^ + [i^ 1 + ^\ 2 ° ~ i^ti <t>\ (50) 

We make use of inequalities (j43|) and (|48|) obtaining that 

||H 2 ^-|^| 2 V|| L1 < (51) 

< c [WWl^WU* + H^^h*] W^h* 
+q|0 1 || L2 [HX^ 1 + ll^lli- -1 ] II^IU* 



< c 
+c 



h<*'- 1 V i \\c\\% + e -MSo-*)/fc|| c ||M2a+l) 



,-(S -y 



|(2<T+1) 

le 1 



S (2a- 1)/4||, 



:|2a +e -2 CT (S -^ 



)/R ||c|ir )2 



e -(S -v)/fe|| c 



I (2(7+1) 

If 1 



< CmaxNIcI 



2<r(2a+l) I, ||2a 1 ||_||(2a+l)„-(S -v)/n 



MWJllci^ 



< Cmax{4 2CT+1)2 ,C +1 }^ (5o ^ )/?l . 
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Hence, 



3 

for any v > and for some C > 0. 

In order to estimate the term (fj,i(c))j^z, which does not contain the vector <j> , 
let us denote, for a fixed j G Z, 

^ = (j) 1 - CjUj = ^ C -rnU m . 

We see that fj,i(c) is given by a finite sum of terms like Cj(uj, <j) 1 \<j) 1 \ 2 ' 7 ) and 

Cj(uj, <^ 1 (^> 1 ) 2cr+1 (</> 1 ) 2 ( cr ~ 1 )); therefore we only have to estimate such a kind of 
term. Indeed, by Lemma|}i. 



^ ICmlllUjUmlU 1 



CfT 



for some v = 0(v). Then it follows that 

< \c m \e-U s °-^-»V H < CMpe-H 8 *-")-*!/* 

m 

< C\\c\\ i xe- < - s °- /)/H < C5 Q e~ {Sa - u)/h 
where we set v + v — > v and ||c||^i < The Lemma is so proved. 

Note that, 

II ||2er+2 in->(j) 1 1 2ct+2 || ||2ct+2 r< 
ll%ll2<7+2 = II TU;u o|| 2( t+2 = IN|| 2( x+2 = °0 

where T is the translation operator defined in Remark [3] Thus, the system (|41j) 
takes the form 

Ac, = A lC ,-/3(rc) J + ( J Dc) 3 + 7 Co|c,f%+7/,(c,0 ± ) , , 

A^ = i^ ± + 7 ii ± |<^ ' { ' 

Dividing by /3, the first equation of (|52[) may be re-written under the following 
form; 



□ 



1 ,^ 



^• = -(Tci, H--(Dc), !- //Icj-l^cy + ^lc. 

where the remainder term ^fj(c, <f>^) is estimated in Lemma H and where D is 
estimated in Remark [5] The approximated system of equations obtained by ne- 
glecting the remainder terms ^(Dc)j and Zfj(c, (j) ) in the Eq. (153")) thus takes the 
form 1|35[) and it is a discrete nonlinear Schrodinger equation equivalent to the so 
called Bose-Hubbard model. 

We are now in position to complete the proof of Theorem [TJ 



J j ' p\—~IJ r 'H^3\ ^3 r nJ3\^i< 



(53) 
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Proof of Theorem [0 Recall that we look for a real- valued solution (j> of the 
stationary equation p4[) . From (|53j) we consider the following mapping 

J-:4xM^4 (c,7)-> 7-(c, 7 ) 

such that 

J"(c,7) = GFj-(c,7)) iez 

and 



F j (c,i) = Ec j + (Tc) j -±(Dc) 



V\cj\ 2a c 



1 



3 - fl/j( c '^ ) i 



where = 0(c, ?i) is the solution of equation (41) given in Theorem[5J This map 
is well defined; indeed we have already seen in Lemma [5] that Uj can be chosen to be 
real- valued, thus for celj, fj(c, </>^) takes real values by construction. This map is 
indeed C 1 in (c, 7) for any a > 0; here, we show that the map f:!gX H 1 — > £^ is C 1 
in c. In order to prove it, we consider, at first, the map c£<j4 fi(c) = (fi,j(c)) ■ 
where fij is defined in the proof of Lemma [9] as 

f 1 , j (c) = {u j ,\4> 1 \ 2 °ct> 1 )-C \c j \ 2 °c j . 

For any c, h G ik we define the linear map from &, to £i 



^i(h) = [{2a + 1) V( Uj , l^ 2 " u k )h k - C \ Cj \ 



By computation we directly have that Z?cfi(c) = A\. For what concern the 
boundedness of A\ we have that the same arguments given in the proof of Lemma 
[9]and the asymptotic estimate of HujUfcHz,! give that 

\\M\c{*^) < C a {l + \\c\\% + ||c|||f)e-^-^/ s . 

In particular Dcfi(c) is continuous in c. Concerning f2(c, <f> ) = (/j,2(c, (p^^j^z, 
which is defined by (|49[) . the same argument applies, and where the map c 1— >• cf) = 
4>(c, h) is smooth by means of the implicit function theorem we have used in the 
proof of Theorem [5] Therefore we have 

.l|2 




\D c f(c,<p)\\ L 



(t 1 



(54) 



From Hyp.2 it turns out that 3 = ^, where C = ||u ||^++ 2 = 0(^ ct/2 ) and 77 
goes to a real value in the multi-scale limit; from these facts and from Lemma [9] 
and Remark [5] wc have that 



= 



and 



(55) 



(56) 



for some < a < So and C > 0, any < v < Sq and any c e ^r(Z) such that 
HclU 1 < <5o fixed. Hence, 



sup 

l|C|Ui<tfo 



Jf(ci(c,fi)) + i(flc) 



- 0(e-«/ 2ft ) =: C h 
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Now, we fix So > 1 and ji% := e a / 2h ^ and we define now the regular mapping 
0:4xR->4, (c,y)^g(c,y) 



such that 



and 



where 



satisfies 



0(c,J/) = {^j(c,J/)}jez 
0,-(c, y) = Ecj + (Tc) . - ?7|c i | 2o 'c j + yg,-(c, ft) 



ffi(c,ft) := — - 
MR 



^/i(c,0(c,«)) + i(5c) i 



sup ||g||^ = 0(1) , asft^O, g(c,ft) = {gj{c,h)} jez - 

\\C\\ e i<8 

By construction at y = then Gj (c, y) coincides with J~j(c, 7). The map C? satisfies 
(?(F, 0) = 0, where F = {Fj}j e % e 4(^) i s a non-trivial real- valued solution of (1351) . 
G(c,y) is C 1 in a neighborhood of (F, 0), with the same reason from (|54|) . for any 
a > and the linearized map 

D c 0(F,O):4->4 

can be written as 

D c 0(F,O)(v) =i+v 

and it is one-to-one and onto (see Hyp. 3). In fact the neighbourhood of the point 
(F, 0) in which C 1 property of Q(c,y) is held is independent of H; let c = F + a 
with ||a||^i < r, and \y\ < S. Then, for any fixed small h > and any e > 0, 



\\D C Q(F + a,y)h- D c G(F,0)h 

< \V\\\{(\F 3 + \F j \^)h j } j \\^ + lylJ-IKlDc/^F + a)^ + ^(Dh)^ 

M>h P P 



< C"|tj|r||h|Ui + \y\— C(l + r 2 + r 2 °)e- a ' 2h \\H^. 
fJ-h 

Here we can take r = C" e/2|rj|, <5 = C~ 1 e/6 (independent of h). 

Therefore, by the Implicit Function Theorem, there exist an h- independent S > 
such that if \y\ < S then there exists a unique solution c(y) in a ^-neighborhood of 
F satisfying Q(c, y) = 0. Since [i% is of order Cn as H goes to zero, y = [i h is in the 
neighborhood for sufficiently small h. Remind that the solution of Q(c, y) = at 
y = [in coincides with the solution of J-(c, 7) = 0. Then we can conclude that there 
exists h* > such that for any h < h* there is a unique solution c £ Bc R (F, 
satisfying J-(c, 7) = 0, for any 7 in a neighborhood of zero. Indeed it follows that 
the radius of the ball (neighborhood) of F is of order Cn since the map y — > c(y) is 
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Then the Theorem follows since 



Y,F jUj \\ H1 < 11^11^ + 



< 



Iff 1 



Efe -Fj) u i 



ff 1 



ff 1 



< II^II^ + Hc-FII^IIuoIIh 1 

< Ce^ So ~ v),h + Ce~ a/2h < Ce~ a/2h 

for any < v < So and for some < a < So, which is exponentially small. 
Theorem [T] is so proved. □ 



3.2. On the validity of Hyp. 3 and proof of Corollary 1. We start with a 
result concerning Hyp. 3 

Lemma 10. Let a > 0. Then Hyp. 3 is satisfied for \rj\ large enough. 

Proof. As a first step we prove that L + is one-to-one by proving that for any 
77* > 1 there exists r\ < with \rj\ > 77* such that 

lli+vll/i > i||v||*, We 4 

Indeed, let F = F(rj) be the real- valued ground state solution of equation (|3"5|) . Let 
G°° be a solution of the following equation 

-\Gf\ 2a Gf = -Gf, (57) 

where G°° € &, ||G°°||#> = 1 for any p e [1, +00]. In fact, we have a family of 
solutions G°° = G°°' J0 for any jo £ Z where G^' 30 — ±5]° ■ Hereafter we fix a given 
value for jo and we denote by G°° the corresponding solution. Recall Remark [TU] 
and we have a solution F = F(rj) of equation ([35)) such that || F(r?) — G°°||^2 < C| | — 1 
and -^j 2 - — > 1 as \rj\ — > 00. Define the following operators 

(L+vJj = -i.( t , J+1 +„ i _ 1 )+(^ -(2a+l)\F j \^)v j 
(ioov), = (l-(2 C 7 + l)|G J -| 2 > i . 
Then, we see, 

IlioovlU, = J2 \ v i - ( 2(J + l )\Gf\ 2 °^\ = E + I - 2av io\ > min {l, |2a|}||v|| £l . 

On the other hand, 
||(L + -L 00 )v||,i < ^— (l^l + ^l +2\vj\) 

+\^r - x | E m + ( 2ct + !) E n^i 2ff - i G ri 2CT iM 
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The first term on the right hand side is bounded by t^tIM^ 1 which tends to zero 
as |ry| goes to oo. The third term is estimated as follows: if 2a > 1, 

^H^r-iGfriK-i < (£\\Ftf° -\GT\ 2 °\ 2 ) 1/2 h\\v 

j 

< c'dlFi^r 1 + HG 00 !!^- 1 )!^ - G°°|H|v||p. 

If < 2cr < 1, 

m\ 2 ° - \Grn\vj\ < ||F - G°l^||v||^ < ||F - G°°||^||v|| £l . 

Therefore, for any rj with I77I > 1 

ll-^+vll^ = \ri\\\L + v\\ e i > \\L + v\\ e i > \\LoovWei - 11(1+ - L oc )v|| <I i 
> (mm{l,\2<j\}-a v -b v )\\v\\ e i 

for some 0<a„(),<l if \r)\ ^> 1. 

By means of the previous result we have that L + is one-to-one since Kcr(i + ) = 
{0}. In order to prove that L + is onto we should require the estimate on the adjoint. 
More precisely, let L* + : £^> — > £^ be the adjoint of the linear and bounded operator 
L + : £g — > ^g. We prove that for any 77* > 1 there exists 77, with \rj\ > 77*, such 
that 

||£+v||*« > ~||v||/=o, VvG^. (58) 

From such an estimate it follows that (see Theorem 11.19 by [5]) L + is onto. 
In order to prove ([58| , first of all we remark that L + acts as 

k 

where a^k = if \j — k\ > 1, 

= OLj,j-i = -1 and oyj = E + 7/(2(7 + l)|Fj| 2<T . 
Since a^fc = ctk,j are real valued then the adjoint L\ formally acts on as 

( L + V ) 7 = J2 a 3,kVk- 



We denote, as before, G°° the solution of equation (j57|) and, similarly, we define 
the operators L + and on ^|P. Then, we see, 

||A»v||*x> =sup|u J - - (2ct + 1)|G°°| 2<t ^| = sup \vj\ + \ -2av jo \ > min{l, |2<r|}||v||/. 
On the other hand, 

|| (Z+ -Loo)v||^ < sup -^(K+il + \vj-i\ +2\vj\) 
iez \V\ 



E-2 

1 



I "'/I 



sup 1^ j + (2(7 + 1) sup \\F 3 | 2CT - |Gf | 2ff \\ Vj I 
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The first term on the right hand side is bounded by A||v||.£°° which tends to zero 
as 1 77 1 goes to 00. In order to estimate the third term we remark that: if 2cx > 1 
\\F.j\ 2a - |G°°| 2<T | < C\Fj - Gf\ [\F 3 \ 2a - 1 + \Gf\ 2a - 1 ] 

2(7—1 1 II f^t 00 II 2a— 11 



< C||F-G°1<- [||F|||r 1 + ||G c 



hence, 

supH^r-IGfriKI < ||F-G°°|U~ [||F||£- 1 + ||G °||&- 1 ] ||v|| 
and if < 2a < 1, 

|2er |/-fOO |2(T 1 1 , | ^ utti /-ioo||2<7 



Bupll^r-IGfHKI^HF-G 

iez 



v 



where ||G°°|U°= = 1, ||F|| £ oo -> 1 and ||F - G°°||^oo -> as |?,| -> +00. Therefore, 
(1581) follows. □ 



Here, we are ready to prove the Corollary. 

Proof of Corollary 1. In order to prove the Corollary we only have to check 
that the linearized map L + : £^ — >■ £^ is, for 77 large enough, one-to-one and onto. 
This fact is proved in Lemma 1101 Hence, the map L + is invertible and we may 
apply the same arguments of the proof of Theorem [T] where, for \rj\ large enough, 
the solution F of equation (|55|) is close to the solution G°° of equation (|57|) . which 
is fully localized on a single lattice cell (see Remarks 1^1 and [TU)) . The Corollary is 
so proved. □ 
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